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Abstract 

Using the supersymmetry technique combined with the transfer matrix ap- 
proach we calculate different physical quantities characterizing localization 
in disordered wires. In particular, we analyze the density-density correlation 
function and its moments and study effects of an external magnetic field H 
on tails of wave functions. At zero and very strong magnetic fields, we obtain 
explicit expressions for all moments and for the entire distribution of the wave 
functions. The crossover between the two limiting cases is more difficult and 
calculations are performed for weak magnetic fields only. We found that the 
far tail of the average density-density correlation function and of its moments 
is strongly influenced by a weak magnetic field and decays twice as slow as 
their main body. Extending Mott's physical picture for the localized states 
we present also a qualitative description of the crossover in the magnetic field. 
From both the analytical calculations and the qualitative description we argue 
that the slower decay of the averaged quantities is a consequence of rare but 
large splashes of wave functions at large distances. The distribution of the 
logarithm of the wave function should also be affected by the magnetic field. 
The splashes and the two-scale decay of the averaged correlation functions 
can be relevant for the conductance of a system of disordered wires connected 
in parallel. 
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I. INTRODUCTION 



Localization of all states in one-dimensional chains and quasi-one-dimensional wires with 
a finite thickness is a very well known phenomenon in the theory of disordered metals. 
Following the first wor kJi where the localization was predicted, a number of analytical 
methods has been developed to treat both the chainslri and wire quantitatively. A 
detailed discussion of the localization in the chains and wires can be found, e.g. in Refs. 

The localization of electron waves in weakly disordered chains and wires occurs due to 
destructive interference and can be destroyed at finite temperatures by inelastic scattering. 
The phenomenon of the electron localization in thick wires is richer than in chains. The 
localization length L c in chains is of the order of the mean free path I. Electrons move 
at distances smaller than the localization length ballistically without being scattered by 
impurities and get localized at distances exceeding L c . In contrast, the localization length 
L c in wires is larger than the mean free path I by a factor proportional to the number of 
channels of the transversal quantization N. At distances larger than I but smaller than L c , 
the electrons diffuse and become localized only at distances larger than L c . 

One more important difference between chains and wires manifests itself if an external 
magnetic field H is applied. In chains the magnetic field cannot influence the electron motion, 
whereas the motion in thick wires is quite sensitive to it. It turns out that the magnetic 
field does not destroy the localization of the wave functions but changes the localization 
length. Remarkably, the localization length L c in thick wires exactly doubles when applying 
a sufficiently strong magnetic field.i'0 

For analytic calculations for thick wires a supermatrix cr-model0lll proved to be a powerful 
tool that allows one to reduce the calculation of kinetic quantities at arbitrary frequencies u 
to solving of a system of differential equations. The localization length L c is obtained from 
an exponential decay of the density-density correlation function. The most interesting limit 
oj — > can be investigated explicitly and one obtains relatively simple formulae for different 
physical quantities proving the localization of all states. 

Calculation of the transmittance of a finite sample with ideal leads is the basis of many 
numerical investigations of the localization properties. The equivalence of different meth- 
ods of calculation of the localization length L c is traditionally guaranteed by the Borland 
conjecturelii stating a complete independence of the localization properties of the boundary 
conditions. 

The limiting cases of a strong magnetic field and no magnetic field, the unitary and 
the orthogonal ensembles, respectively, are relatively well studied in quasi ID, whereas the 
crossover between these two limits has not been well understood. It was assumed that the 
localization length changed smoothlyBEl between L c and 2L C and one had to calculate a 
curve connecting L c and 2L C when increasing the magnetic field from to oo. 

At the same time, investigation of the crossover can have quite interesting applications 
now because a systematic experimental study of the localization has been performed recently. 
In the experiment, a large number of submicrometer-wide wires prepared from Si 5- 
doped GaAs were connected in parallel to insure statistical averaging. This enabled accurate 
measurements of the conductivity as a function of the temperature T and the magnetic field 
H applied perpendicular to the wires. The regimes of the weak and strong localization were 
observed and the experimental data in the regime of strong localization were used to extract 
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the dependence of the localization length L c on T and H. The localization length L c was 
shown to double in a strong magnetic field as compared with the limit of zero magnetic 
field H, which is in a complete agreement with the theory.!'! As the measurements are done 
at an arbitrary magnetic field, study of the crossover between the orthogonal and unitary 
ensembles may become verv important. 

In a recent publication,^ we discovered that already an arbitrarily weak magnetic field 
drastically changes the tails of the averaged amplitudes of the wave functions. At suffi- 
ciently large distances depending on the magnetic field, the exponential decay with the 
localization length L c changes to a decay with the length 2L C . This result is of interest not 
only from the theoretical point of view but it might also be important for the transport 
measurements,0Jill since the localization length L c is relevant for the hopping conductivity 
at low temperature.00 

However, a numerical study of the averaged logarithm of the transmittance performed for 
sufficiently long wires0 did not manifest any two-scale behavior showed a smooth variation of 
the localization length between L c and 2L C . No sharp change in the far tail of the logarithm 
of the transmittance was observed in a weak magnetic field. As a possible explanation of 
the contradiction several reasons were suggested. As one of the possibilities, the difference 
between the logarithm of the averaged wave functions, the quantity we have considered^ 
and the averaged logarithm of the wave functions was mentioned.0 

In this paper, we give details of the calculations presented in our short publication!!! and 
discuss the question in what situations our effect of the changing of the tails of the wave 
functions by a weak magnetic field can be identified. To understand better the shape of 
the wave functions we calculate not only the density-density correlation function but also 
its moments. The changes in the tails of the wave functions are seen in all moments of 
the density-density correlation functions and all of them are proportional at large distances 
x to exp(— |x|/4L cu ), where L cn is the localization length for the unitary ensemble. The 
independence of the exponential decay of the asymptotics of the number of the moment 
signals that they are formed by strong rare splashes of the wave functions. 

The situation with the averaged logarithm of the wave function is more delicate. The 
large splashes contribute also to the averaged logarithm of the wave functions at finite dis- 
tances. As a result, some broad distribution of logarithms occurs. Unfortunately, we are 
not able to predict how this distribution function develops in the limit of infinitely long dis- 
tances. However, our results provide us reasons enough to suggest that the logarithmically 
normal distribution implied within the Borland conjecture is not universally applicable for 
length of the sample x ^> L c for the description of the crossover region. Existing numerical 
studiesil'ii'0 might be interpreted as the head of some logarithmically normal distribution 
changes smoothly with H. Our results do not exclude a smooth behavior of the average log- 
arithm of the wave function. We show, however, that another length scale, xh ~ L c ln(l/H) 
enters the theory of localization, up to which the distribution should differ from the stan- 
dard logarithmically normal one. Thus, the whole distribution function of the logarithm is 
of interest. 

For calculations we use the supersymmetry technique that makes possible consideration 
of disordered wires with any kind of boundary conditions. Inparticular, one can study the 
transmittance of the disordered wires with metallic leadsffl An external magnetic field 
is naturally accounted for within supersymmetry formalism allowing study of the crossover 
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between the unitary and orthogonal ensembles. 

The localization of the wave functions discussed in the present paper is extracted from 
the moments of the density-density correlation function 



pW ( X - x >) = (T \i>a (x)\ 2n (x')\ 2n 6(e - e a )\ (1) 

where ip a (x) and e a are the eigenfunction and eigenenergy of a state a, x and x' are the 
coordinates along the wire, and the angle brackets stand for the averaging over the disorder. 
This quantity was previously studied for a closed finite multichannel wire,^3 when the wave 
functions were taken at the ends of the wire. Generalizing the expressions for arbitrary n, 
the complete distribution was also obtained.0 

Naturally, the very presence of hard boundaries or ideal leads inevitably changes localized 
states near the ends of the disordered sample. These states should differ from those localized 
in the bulk. Such a difference has been found indeed for the density of states of a single 
chain. H 

The localization properties extracted from the transmittance of an open system^lS and 
from the density-density correlator of a closed one manifest a generic universality: All 
moments of the both quantities decay with the same rate, exp(— \x — x'|/4L c ), and the 
logarithm of the both quantities is normally distributed. The fact that the decay rate is the 
same for all moments shows that an important contribution comes from rare strong splashes 
of the wave functions. Due to these splashes the moments decay slower with the length 
4L C and not with the localization length L c . At the same time, the logarithm of physical 
quantities like the transmittance or the product of the amplitudes of the wave functions 
at different points is a self-averaging quantity and approaches the value —\x — x'\/L c for 
| a; — x'\ — > oo. This quantity characterizes typical wave functions. The question what 
happens with the wave functions if a weak magnetic field is switched on has not been 
discussed. 

In the subsequent chapters we present details of our analysis. We use an approach based 
on the Lebedev-Kontorovich transformation to find the exact distribution of any quantity 
of interest in quasi- ID wires for the orthogonal and unitary ensemble. This allows us to 
calculate the correlator, Eq. ([]]), and the Landauer-like conductivity. Exact expressions 
valid at arbitrary distances are obtained for all moments for the pure orthogonal and unitary 
ensemble. Using the formulae for the moments we derive the entire distribution function. At 
a weak magnetic field, we show that a peculiar two-scale behavior characterizes any moment 
of the wave functions, Eq. (|]). 

The paper is organized as follows: In Sec. II we introduce the quantities to be studied. 
A transfer-matrix formalism used for the calculations is described in Sec. III. Reduction to 
an effective Schrodinger equation is carried out in Sec. IV. Then, in Sec. V we calculate the 
moments and the distribution functions in the limiting cases of the orthogonal and unitary 
ensembles. In Sec. VI we present results for the wire in a weak magnetic field. Using a 
standard qualitative picture for localized statesillll we discuss the physics involved in Sec. 
VII. Our arguments allow us to obtain without any calculation the characteristic length 
scales of the problem and reach a better understanding of the localization in the magnetic 
field. We discuss the results and make conclusions in Sec. VIII. 
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II. CORRELATION FUNCTIONS 



Localization in disordered metals can be well characterized by different correlation and 
distribution functions. For an infinite sample one can consider the density-density correlation 
function and its moments. For a finite wire connected to metallic leads the transmittance 
and its moments are of interest. 

These quantities can be efficiently studied using the supersymmetry method. Follow- 
ing this approach one should reduce calculation of the correlation functions of interest to 
computation of correlation functions of a supermatrix cr-model. In the standard formula- 
tion the cr-model contains 8x8 supermatrices Q. Fluctuations of the supermatrices Q are 
strongly influenced by a magnetic field and one half of the "degrees of freedom" is sup- 
pressed if the magnetic field is sufficiently strong. This corresponds to reducing the size of 
the supermatrices Q to 4 x 4. Keeping this size of the supermatrices is sufficient to calculate 
such physical quantities as conductivity or density-density correlation function. 

In order to calculate higher moments of these quantities or distribution functions one has 
to increase the size of the supermatrices, which would make the theory very complicated. 
Fortunately, calculation of moments of certain quantities does not demand increasing the 
size of the supermatrices, which allows to compute them explicitly. 

In all such cases the free energy functional F [Q] describing fluctuations of the 8x8 
supermatrices Q can be written as 



F[Q] = —Sir J D ( V r Q(r) - ^A[Q(r), r 3 ] + 2iuAQ(r) 



dr. (2) 



The cr-model, Eq. appears after averaging over disorder and we are to carry out com- 
putations with the regular model. In Eq. (fj) D is the classical diffusion coefficient, v is 
the density of states, and uj is the external frequency. The supermatrix Q obeys the con- 
straint Q 2 = 1, A is the vector potential corresponding to an external magnetic field H, 
and the standard notations for the supertrace Str and matrices A, r 3 are used. El Any quan- 
tity of interest in the present formulation is expressed as a functional integral over Q(r) 
with the weight exp(—F[Q}) and a combination of elements of the supermatrix Q in the 
pre-exp onent ial . 

We study properties of wave functions calculating their moments written in Eq. ([!]). As 
the first step, we express the moments of the wave functions in terms of retarded G R and 
advanced G A Green functions that can be written as 

G^M = ^ a{X) f AX '\ (3) 

6 E a 

a 

where ip a and e a are eigenfunctions and eigenenergies of the electron states in the disor- 
dered system. Using the spectral expansion, Eq. @, assuming that all states are localized, 
such that the spectrum is discrete, and considering the most divergent terms at vanishing 
frequency uj — > one can come to the following relation 

[Gf (x,x')G* +e (x',x)] n = 2m l 2n ~^ ( "J; ) 7 1 E \Mx)\ 2n \Mx)\ 2n 5(e - e a ) (4) 

((n - 1)!) w 2 ™- 1 
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Although Eq. (|j) contains an arbitrary power of the Green functions we can express the 
LHS in terms of 8-component supervectors. This possibility is due to the fact that the Green 
functions are taken at two different points only. The next step is the averaging over disorder 
and the decoupling the effective "interaction" by integration over the supermatrices Q. A 
functional integral over Q is simplified in the saddle-point approximation and we come after 
standard manipulations^ to the following expression 

([G*(x,x f )G* + M^)] n ) = (™) 2n ([Qll(x)} n [Qll(x')} n ) F . (5) 

In the LHS and RHS of Eq. ([5]), averaging over the impurity potential and averaging over 
the free energy, Eq. (|]), is implied, respectively. The superscripts and subscripts of the 
supermatrices Q stand for certain matrix elements. 

In Section V, exact expressions for the moments of the density-density correlation func- 
tion, Eq. ([I]), as well as its entire distribution function 

V^t) = (j^6(t- (x)\ 2 \ij a Or')! 2 ) 5 (e - e a )J . (6) 

will be obtained. 

Hereinafter we use dimensionless expression for the wave functions ipSL c — > if), where S 
is the cross-section of the wire and L c is the localization length. To calculate the distribution 
function V^(t), Eq. (||), we introduce an auxiliary correlation function 

V a , b (t u t 2 ) = (5 (h - aQll( Xl )) S (t 2 - bQH(x 2 ))) F , (7) 

where a and b are some parameters. Once the correlator Eq. (|7|) is known, any other quantity 
of interest can be found by integrations of the type V a ,b(a) da using a proper choice of 
a, b(a). As an example, we will calculate the function and the distribution 

V Q = (5 {t+( i u;L c uS) 2 Qlt(x)Qll(x'))) F , (8) 

characterizing the Landauer-type conductivity.^ Hereinafter, all the correlators are studied 
in the limit to — > 0. 

To express the distribution Eq. (^J) in terms of a functional integral over the supermatrices 
Q, we relate the numerical coefficients of all moments of the distribution functions and 
Vq, using the identity 



(2n-2)!(-l) n d 
((n-1)!) 2 dp 



(3=1 JO 



[ dp^r-^i-p))"- 1 . (9) 

Jo 



Thus, we reduce the entire distribution to Vq 



Eq. fllPp can be proven taking (t) from Eq. (0) and expanding both sides in ip a and 
Q respectively. Then, one should use Eqs. (f|, |5]) and this gives finally Eq. ( JTOj ) with the 
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functional F from Eq. (0). Eq. ( |10D allows us to calculate the distribution function of the 
amplitudes of the wave functions at two different points of an infinite sample. 

Although we calculate here the distribution functions of the quantities written at two 
different points x and x', our procedure can be extended to many-point functions. 

In quasi-lD case, one can calculate the functional integral, Eq. (|7|), exactly using the 
transfer-matrix technique. Within this method the functional integration is reduced to solv- 
ing an effective "Schrodinger equation" in the space of variables parametrizing Q-matrices.i 
Solving the "Schrodinger equation" and calculating certain "matrix elements" that are def- 
inite integrals over Q one can obtain any physical quantity of interest. For instance, the 
moments corresponding to the distribution function V a ,b, Eq. (0), read 

T mn = (%t™) VQ = a n b m J *(Q) (Qlt) n T(x,x';Q,Q') (Q%) m *(Q') dQdQ', (11) 

where the function \l/ (Q) is the partition functions of the parts of the wire to the left from 
the point x and to the right from the point x', and the function T represents the partition 
functions of the segment between x and x'. 

The meaning of the functions \I/ and T becomes clear after the reduction of the functional 
integrals to the effective Schrodinger equationBlll'ill The function \1/ (Q) appears to be the 
ground state of the effective Hamiltonian Tiq acting in the space of the elements of the Q- 
matrix whereas the function V is the Green function in this space. They satisfy the following 
equations 

n Q ^ (Q) = (12) 

+ hq ) r (a; ' x ' ] Q1 Ql) dx = 6 ( x - x ') 

The fact that the energy of the ground state is zero is a consequence of the supersymmetry 
of the initial a- model, Eq. (|2]). 

The explicit form of the Hamiltonian 7{q depends on the parametrization of the Q-matrix 
and its choice depends on the quantities calculated and the physics involved. 

The supermatrices Q consist of compact and non-compact blocks containing circular and 
hyperbolic functions, respect ively.0 In the cases of zero and strong magnetic fields only the 
non-compact sector survives in the limit u — > 0, which greatly simplifies the whole analysis. 

If a non-compact variable Ai is relevant only, the function \1/ takes in the limit of small 
cu the well-known formM 

^ = 2z 1/2 K 1 (2z 1/2 ), z = (13) 

valid both for chains and thick wires. Using the modified Bessel functions Ki P one can write 
in this limit the solution for the function T as 

o poo 

T(x, x'; z, z') = 4 / e - |x - x ' l(1+p2)/4ic z 1/2 KJ2z^ 2 ) z' l/2 K ip (2z' 1/2 ) sinh(7rp)pdp . (14) 

Jo 

where L c is the localization length. Eq. (IT3) will be proven in Sec. V. Actually, it follows 
from the orthonormality of the modified Bessel functions Ki P (see Eq. (|42|)) and the fact that 
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the functions 2z 1 / 2 K ip (2z 1 / 2 ) are eigenf unctions of the "heat" equation with the Laplacian 
in the space of the matrices Q. 

The expressions for \1/ and T, Eqs. (13, [14]), are sufficient to calculate all the moments, 
Eq. (|TTD , for the orthogonal and unitary ensembles. This way is more straightforward and 
more economical, as compared to the usual of representing the moments as solutions 

of differential equations obtained within the transfer matrix technique. At the same time, 
the conventional approach is more general and is applicable also in the crossover regime. 

It is worth mentioning that our analysis of the localization properties in infinite wires is 
somewhat easier than that for conductanceBU and wave functions!! in finite wires. This is 
because only large values of the non-compact variables Ai are important for the calculations 
in our case. This is this simplification that allows us to calculate all moments as well as the 
entire distribution functions for the 'pure' (orthogonal and unitary) ensembles (see Sec. V). 



III. TRANSFER-MATRIX TECHNIQUE 

At not very high frequencies, the zero transversal space harmonics of the supermatrix 
Q in Eq. (fj) gives the main contribution. In order to neglect the non-zero harmonics one 
should choose the London gauge for the vector potential A. Writing A = (Hy, 0, 0) , where 
H is the external magnetic field perpendicular to the wire we represent the free energy, 
Eq. (0), for the quasi- ID sample in the form 

F[Q] = ^Str J 

where L cn = is the localization length for the unitary ensemble. The crossover 

parameter X depends on the geometry of the sample, X = 27r0/0 o , 0o = hc/e is the 
flux quantum, and = HL cn {y 2 )^ 2 is the magnetic flux through the area limited by the 
localization length. The brackets imply averaging across the wire over the coordinate y 
which is chosen to be perpendicular both to the direction of the wire and the magnetic field. 
This averaging gives (y 2 ) 1 ^ 2 = dj y/l2, where d is the width of the wire, for "flat" wires made 
on the basis of a 2D gas.00 For wires with a circular cross- sect ion, (y 2 ) 1 ^ 2 — d/A, where d 
is the diameter. 

To determine the functions \1/ and T in Eq. (|TTJ) , we discretize the wire by introducing 
sites on which the Q-supermatrix varies. The function \1/ (Q) is the partition function of a 
semi-infinite sample and therefore it does not depend on the coordinate x. This allows us 
to write immediately the a recurrence equation. Relating the function ^{Q) on one site and 
its value \1/ (Q f ) on the neighboring one we obtain 

*(Q) = J N{Q,Q')Z {Q')*(Q')dQ' . (16) 

where the function N (Q, Q') describing the coupling between the sites 

N(Q, Q') = exp ^|^|Str QQ'^j , (17) 

originates from the kinetic term (VQ) 2 in the free energy, Eq. (p~5|) . The function Z (Q), 



(V a Q) 2 + 



X 2 



IQL 2 



■lQ,Tsf 



2iuAQ(x) 
D 



dx . 



(15) 
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MQ) = exp (^g^Str AQ + ^Str [Q, r 3 ] 2 ) (18) 



describes the remaining terms in the free energy functional F, Eq. (|15|), and should be inte- 
grated on each site. In order to reduce the integral Eq. flUf ) to the differential "Schrodinger 
equation" for the functions \l/ (and similarly for T), one should take the continuous limit 
Al — > 0, where Al is the length of the separation between the sites. In this continuous limit, 
Al drops out from final expressions and therefore can be chosen to be arbitrarily short. 

Since the kernel T enters equations for each moment only in a combination with 
(Q 2l ) n ^(Q), it is advantageous to introduce a matrix function Pjf 1 ' such that 

Pt\Q) = J2 ex v^ x - x 'W I mQ,Q')(Q' 21 mQ')dQ', (19) 

x J 

where T(k; Q, Q') is the Fourier transform of T(x — x'; Q, Q'). The matrix function P k n) (Q) 
on a site is related to its value P k n) (Q') on the neighboring site as 

P { k n \Q) -exp(-ikAl) J N(Q,Q')Z (Q')Pl n) (Q')dQ' = (Q 21 ) n V(Q) . (20) 

The moments of the distribution Vq, expressed in terms of the functions ^ and P k , are given 
by 

T mn = J HQWlYiPlt + P ( -LW , (21) 

Solving Eqs. ( |19| - p0] ) and calculating the integral in Eq. ( |21| ) one can find, in principle, all 
the moments. However, further simplifications are necessary to obtain the results explicitly. 

A spectral expansioni'Elii is convenient for the analysis of the function P^- So, we 
expand the functions P^ in eigenfunctions 4>s(Q)i introduced as 

J N(Q, Q')Z (Q')MQ'W = £MQ) ■ (22) 

The matrix functions <ps{Q) have the same structure as Q 21 . Their orthogonality and the 
normalization properties can be written in the form 

J [<i>t>(Q)U[MQ)U*Q = s (£ - E ') ■ (23) 

Expanding P^ in the complete set of the functions 0£, Fourier transforming it to the 
coordinate representation, and substituting the result into Eq. ([0]), we arrive at 

T mn = -^exp(-^|x-x'|) J(Qll) m [MQ)h^(Q)dQ J(Qll) n [MQ)h^(Q)dQ. (24) 



Eq. ( p^ ) contains the sum over all eigenstates that should be found from Eq. (22[) and 
establishes a link between the eigenvalue problem in the space of parameters of the Q- 
matrix and the localization properties. The lowest non-zero eigenenergy corresponding to 
the first excited state determines the exponential decay of the correlation functions at large 
distances. 
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IV. EFFECTIVE HAMILTONIANS 



Eqs. flTB], E4T) form a closed system completely solving the problem. To perform 



explicit calculations a parametrization of the Q-matrices should be chosen. For theDrpblem 
involved one can use either the "standard" @Jl3 or the "magnetic" parametrizationSlI. The 
former parametrization has been used, e.g., for calculation of the density-density correlation 
function in disordered wires for the orthogonal and unitary ensembles, while the latter one 
helped to obtain explicit formulae for the crossover between the ensembles in the zero- 
dimensional situation. 

For the problem of the localization in wires the standard parametrization allows to 
calculate any quantity of interest for the orthogonal and unitary ensembles and it will be 
used in Sec. V for calculations of all moments as well as of the distribution function of the 
density-density correlator. At the same time, the standard parametrization is of little help 
in the crossover regime between the ensembles and we will use in this case the magnetic 
parametrization. 

The basic equations and their solutions for the orthogonal and unitary ensembles can be 
found in Refs.B@ and we do not write them here. Instead, let us present explicit differential 



equations that can be obtained within the magnetic parametrization from Eqs. ([16], |22|, ^4|) . 

The Q-matrix in this parametrization is constructed from two sets of variables, cooperon 
ones being cut by the magnetic field, and diffuson ones insensitive to H. An explicit form of 
the supermatrices Q in these two parametrizations can be found elsewhere. §The supermatrix 
Q in this parametrization is represented as 

Q = V d Q c V d , (25) 

where the supermatrix Q c , Q 2 C = 1, contains cooperon "degrees of freedom" and V d , V d V d = 1, 
diffuson ones. To make the reading easier, we present here the essential structure of the 
block Q 12 in the magnetic parametrization, keeping only those terms that contribute to the 
correlators considered 

Q 12 = u d i{l - \ d ) 1,2 [k + 2(r, c r,* c - k c k* c )(\ 1c - X c )]v d , (26) 

where X C)d = diag(cos# Cj d, cosh#i C) i^), v d and u d are the standard 4x4 unitary matrices in 
the V d (diffuson) block, and T] c , tj*, k c , k* are the Grassmann variables from the cooperon 
block. 

An advantageous feature of the magnetic parametrization is that the term containing 
the magnetic field H in the free energy has a simple form 

Str[Q,r 3 ] 2 = 16(A 2 C -A 2 ). (27) 

In contrast, the term with the frequency to is quite complex since it contains the Grassmann 
variables 

Str(AQ) = 4[A ld A lc - \ d \ c + 2( VcV * c - k c k* c ) {X lc - A c ) (A w - X d )) . (28) 

The functions \l/ and (pg correspond to states having respectively the "angular momen- 
tum" zero and one with respect to the unitary rotations V d . The form of the functions </>£ 
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is determined by the RHS of Eq. fl2~0|) and, as a result, it should be searched for in the class 
of functions with the same structure as the matrix Q 12 

^O^^al-Va)*- (29) 

where Vd and Ud are unitary matrices containing Grassmann variables only The different 
forms of the functions \I> and (j)g lead to different forms of the effective Hamiltonians de- 
termining the functions \I/ and (j)g. The functions / and f\ do not contain the diffuson 
Grassmann variables but may contain all other variables. 

The Schrodinger equations corresponding to Eqs. (|1(J) and fl22| ) are obtained considering 
slow variations of Q on neighboring sites. The matrix Q' is expanded near Q and integrated 
over small variations of the matrix elements. The expansion up to quadratic variations is 
sufficient and therefore, the quadratic form obtained is identical to the one describing the 
square of the elementary length in the space of the supermatrices This elementary length 
determines the Berezinian corresponding to a chosen parametrization. Correspondingly, 
the "Laplacian" entering the effective Schrodinger equation contains the Berezinian of the 
transformation. The Berezinian of the magnetic parametrization readJ^E 

1 AX 2 C 

J = J c JdJcdi Jc,d — ofi 2 / , > To" , Jed = 77 j TTo") (30) 



The parts and of the Berezinian Eq. ( [30]) originate from the diffuson Grassmann 
variables rjd, Kd, whereas J c comes from the cooperon ones r] c , k c . The diffuson Grassmann 
variables do not enter Eqs. fl2"7|) and (|28|). In some sense, the diffuson variables decouple 
from other variables and this simplifies calculations. 

The situation with the cooperon variables is more difficult. Since the free energy, Eqs. ( TT5| , 



28|) , contains explicitly the cooperon Grassmann variables r\ c and k c , so do the functions \l/ 
and <p£. In order to write equations for these functions explicitly one should expand them 
in powers of the Grassmann variables. Proceeding in this way one obtains a system of 
differential equations for the functions entering the expansion. Another complication occurs 
due to a non-trivial angular dependence on Q c . This matrix is surrounded by the diffuson 
modes and couples to them. As a result, one comes to extremely cumbersome equations 
that can hardly be solved analytically for arbitrary X. The solution can be found for X ^> 1 
but this limit corresponds to the well studied unitary ensemble. 

Fortunately, an important information about correlation functions can be extracted in 
the limit I < 1. To understand why this limit helps to study the problem we remind 
the reader that the possibility to obtain closed expressions for correlation functions in the 
orthogonal and unitary ensembles is due to the fact that the main contribution to integrals 
comes in the limit of small frequencies u from large values of the non-compact variables (in 
the unitary ensemble Aid ~ V^)- For large values of X\d the partial differential equations 
for the functions \1/ and (f>s are sufficiently simple and can be solved. In the crossover regime 
when X is finite, the equations are still very complicated even in the limit uj —>■ 0. However, 
in the limit X 1, some quantities of interest are determined by large Ai c ~ 1/X ^> 1. This 
leads to an additional simplification of the equations and the possibility to make estimates. 
At the same time, not any correlation function can be calculated in this way, which limits 
the applicability of the trick. 
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(31) 



In the continuous limit, Eq. ( jig) can be reduced to the form 

Hx^ = 

Due to the presence of the Grassmann variables in the free energy functional, Eqs. (|T5| , 
8]), the solution \1/ must contain them too, which is in contrast to the solutions for the 
orthogonal and unitary ensembles. The function \1/ can be represented as 

* = (Aid, Ai c ) + #i (Ai c , A i(i ) (VcVc ~ ^K) + ^2 (Ai c , Ai d ) ^*k c k* (32) 
Then, one obtains a system of partial differential equations for the functions \l/o; ^i, ^2 

J (Ai c - A c ) 2 



nf + ) *i 



AicA c 



(Aic — A c ) 



(33) 



where 



(nf + uxA ^ 2 - 2ux 2 ^ 1 = 



n 



(0) 
X 



^ -^|1-A-|^- ^ yy^^JcJcdlA, 2 -!!— + X 2 (A? C -A 2 ) 



Z -. J d d\ l 



i=d,ld 



d\i \ J c Jcd d\i 

i=c,lc 



(34) 



and 



Q = 2n 2 (uS) 2 Duj 



xi — AidAi c — AdA c , x 2 — 2 (Aic — A c ) (A^ — A d ) 



(35) 



In the limit Aid, Ai c ^> 1, the difference between equations for the functions \& and is 
negligible because it comes from the averages of the type (ArjiArj*) ~ 1/Aii, (A^Ak*) ~ 
I/Ah, i = c, d. In this limit, Eqs. (|33| , 34] , |35|) take a simpler form 



+ uxi ) ^1 + ux 2 ^ = 



(36) 



(ft?? + ^2 - 2cDX2^1 = 



with 
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- - Xl^- + 2A lc |- + X 2 A? C (37) 



xi = x 2 /2 = XiAid 

The same equations can be written in the limit Aid, Xu ^> 1 for the functions </>£. 



It is clear that Eqs. ( |32| , |35| ) cannot be solved analytically for an arbitrary X and the 
only hope is to solve Eqs. (|36| , [371) . Of course, one is restricted then by calculation of those 
quantities for which large A^, Ai c give the main contribution. 

In the opposite limiting case X —>■ 00 but u — > all the cooperon variables are frozen 
and one comes to the effective Hamiltonian 

d 2 2iuL 2 „ , . , 

n = - x ^-^r x ' (38) 

After a proper change of the localization length the same Hamiltonian can be written in 
the orthogonal case using the standard parametrization. For completeness let us remind that 
the localization lengths for the orthogonal and unitary ensembles, L co and L cn respectively, 
are given by the following expressions!@ 

L cu = lixvSD , L co = ixvSD , (39) 

In these cases the function ^ does not contain Grassmann variables. The equations to 
be solved take the form 

HV = 0, H<j)s = L c S<J )£ (40) 

where the localization length L c equals either L co or L cn depending on the ensemble. 

In the next Section, solving Eqs. ( {4*0]) we will calculate all moments and the entire 
distribution function of the density- density correlations. This will help us to come to certain 
conclusions about properties of the electron wave functions in the orthogonal and unitary 
ensembles. 



V. CALCULATION OF MOMENTS AND DISTRIBUTIONS FOR THE 
ORTHOGONAL AND UNITARY ENSEMBLES. 

Calculation of the moments and of the distribution functions is important for the under- 
standing of the structure of the wave functions and properties of different averages. Higher 
moments and the distribution function for infinitely long wires have not been considered yet 
although the results are known for finite wires£3 So, calculation of the moments and of the 
distribution function can be interesting on its own but, what is more important, this infor- 
mation will also help us to make certain conclusions about properties of the wave functions 
in a weak magnetic field where the possibility of exact calculations is more limited. This will 
help us to extend the results of Ref.@ where only the averaged density-density correlation 
function was calculated. 

Solutions \1/ and <p£ of Eqs. (|3^, |40|) are well known (see e.g.EEB) and can be written as 
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= 2^/2^(2^/2) ) £ = 2z 1 /2^ v ( 2 ^i/2) ; f = (i + p 2 )/4L c , (41) 

where z = 2i\]ioL 2 c j D . Comparing Eqs. (Q) and (|TT|), we can understand that the prop- 
agator has the form of Eq. fll4|) . The numerical coefficient in Eq. ( |i~4"|) is found using the 
Lebedev-Kontorovich transformation^!! applicable for an arbitrary function (p(y) 

m = JT v(v)Kip(y)f , <p{v) = ^ JT v{p)k»(v)p^{*pW (42) 

Now we can calculate all the moments of the two-point correlator, Eq. (0). Substituting 
Eqs. ( f4"T|) into Eq. ([24]), and calculating integrals over all elements of the supermatrices Q 
except Ai, yields 

or) /-00 

Tmn = {t n 1 t%) VQ = — dppsmh('Kp)e-n p2+1 )M n (p)M m (p) . (43) 



where 



M n = y dta f|) rfK^tjKipih, 



and m = |x — x'\/L c . 

The structure of Eq. (^) is not complicated and the integration over the variable t\ can 
be performed exactly^. As a result, we come to general expressions for all moments of the 
density-density correlations p^o (u) introduced in Eq. ([[]) 

(B) 21/ (2n-2)! /"» , (p2+1 



7T 



(("-1)1 



,2 



/*0O 

/ dppsinh(7rp)e-^ 2+1 ) [P (n) (p)] 2 , (44) 



n 2 f 1 ip\ / 1 ip\ ( 1 «p\ / 1 ip 



n 



(2n-l)! V 2 2/ V 2 2/ ^ 2 2/ V 2 2 J 

For n = 1, Eq. ( fill) reduces to the expression obtained by Gogolin0 for a single chain. Eq. 
(fH|) is applicable at any x including the values smaller and of the order of L c for both the 
unitary and orthogonal ensembles. Comparing the moments, Eqs. ([44]) , with those obtained 
for finite wirescZl, we see that their structure is similar, although they are not completely 
identical. 

At large distances, \x — x'\ ^> L c , essential contribution to the integral over p in Eq. ( fPj) 
comes from small p <C 1 and the integral can be easily calculated. We see that the moments 
of all quantities decay similarly as in infinite wires, Eqs. ([Hp, (f4"3"D, finite closed wiresi^, as 
well as in open chainsHS 

p£»W = C„(^) 3/2 exp(-M) . (45) 

where C n is a constant depending on n. The exponential decay of all moments seen from 
Eq. ( fl~5D manifests the localization of the wave functions. A very important feature of the 
moments p£?(x), Eq. (f4~5), is that the dependence on the coordinate is the same for all n. 
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This does not allow to interpret (x), Eq. (|45|) , as functions characterizing the typical 
shape of the wave functions. 

The exact representation of the moments of the wave functions, Eq. fl43|) , allows to obtain 
the entire distribution function of the density-density correlations. Of the crucial importance 
is the fact that Eq. is valid for all n including the limit n — > oo. It is clear, if the nth 
moment of a quantity can be represented in a form of an integral of a function multiplied 
by the nth power of the variable of integration, then this function is the distribution of this 
quantity. In the case under consideration, Eq. (^) contains, besides the necessary power 
of the variable t, the factor n 2 , originating from the expansion of the element {Ql 2 ) n in the 
Grassmann variables. This does not make the calculation of the distribution function more 
difficult because the presence of this factor leads merely to additional derivatives, and we 
obtain for the distribution function V(t 1 ,t 2 ), Eq. (|7|) 

V{t l ,t 2 ) = — / dppsm\i(7rp)e-n p2+1 )U p (t 1 )U p (t 2 ) , (46) 
n Jo 

where 

H,(t) = ([^ (2t 1 / 2 )K ip (2t 1 / 2 )] ' t y . (47) 

and ' stands for the derivative over t. The two-point distribution function V(ti,t 2 ) can 
be used to obtain any other correlator of interest. For instance, the distribution of the 
Landauer-type conductivity Pq(£), Eq. (§) can be represented as the convolution 

V Q = y dr (6{t - tujuLcSQ^MMt - iuuL c SQ%{x 2 ))) F . (48) 

Then, we obtain for the function Vq (t) 

V Q (t) = - dppsinh(7rp)e-f(" 2+1 ) / — n„(r) U p (t/r) , (49) 
n Jo Jo T 

Eq. (f49D is also exact and applicable at any distances. It is interesting to find this 
function near its maximum and this can be done easily in the limit \x — x'\ ^> L c . At very 
large separation between x and x', typical values of t are exponentially small in distributions 
Vq and V^, Eqs. (§, 0). A dominant contribution to the integral, Eq. (|49D, comes from 
the region r ~ t < 1. Eq. ( |47|) allows us to write the asymptotics of the function n p (r) at 
t < 1 and t > 1 

n p (r) ~ Re [T(tp)/4T 3/2+if,/2 ] , r < 1 



n p (r) ~ vrexp (-4r 1/2 ) /r 1/2 , r > 1 (50) 



which ensures convergence of the integral over r at r — > in Eq. ([49D. Performing integration 
in the resulting expression 



_ . , . (°° , . w x T(ip)r(2 + ip) / it ln 2 t m. ilnt N2 
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in the saddle point approximation, we obtain 



^ g (-ln*)~^i^^(-^(« + ^*) 2 ) • (51) 

Considering the distribution V^, Eq. (flQ|) , under the same assumptions it is not difficult 
to check that in the limit \x — x'\ ^> L c it reduces to Eq. 

The saddle-point approximation used for the derivation of Eq. fl5l|) is applicable at large 
distances \x — x'\ ^> L c when u is large. Eq. (|5l|) shows that in this limit the distribution 
functions are log-normal. The result about the log-normal distribution of physical quan- 
tities like conductance agrees with the corresponding result obtained for chainsilS and is 
universal. It shows that the localization length of the wave functions or in, other words, the 
logarithm of the conductance is a self-averaging quantity. Considering very long samples one 
has to obtain the localization length L c and fluctuations of this quantity can be neglected. 

At the same time, the log-normal distribution, Eq. ([51]), is applicable for not very large 
values of t only. The far tail of the distribution function decays more slowly. The 

moments of the density-density correlations Poo, Eq. cannot be obtained from Eq. 



(pl|) because the main contribution comes from the far tails of the distribution function. All 
these peculiarities are discussed in details in RefsJiH. 

What is the physical meaning of the fact that all the moments Eq. (p[5|) decay in the 
same way with a larger length 4L C ? This can be understood if we imagine that there are big 
splashes of the wave functions even far from the localization center. Although, they can give 
a considerable contribution to the shape of the wave function, their probability decreases 
exponentially with |x| (is proportional to exp (— |x|/4L c )). Then, such splashes would give 
a similar contribution to all the moments. As concerns typical values of the squares of wave 
functions, they decay faster with the localization length L c and give a smaller contribution 
to the moments. 

The situation is opposite if one calculates the average logarithm of the wave functions. 
For the logarithm, the splashes do not give a specially large contribution, being at the same 
time very rare. So, their contribution to the average logarithm can be neglected. In other 
words, the average logarithm of the correlations of wave functions characterizes typical wave 
functions, whereas the moments describe the rare splashes of the wave functions. This 
understanding will be important for interpretation of the results in a weak magnetic fields 
obtained in the next Section. 



VI. EFFECT OF MAGNETIC FIELD ON THE TAILS OF WAVE FUNCTIONS 

This section is central in our paper. We try now describe electron wave functions in 
disordered wires in a weak magnetic field, such that J < 1. General equations describing 
this situation are very difficult even when they are written for the "ground state" function 

Eqs. (|3ll - |35D . At the same time, any attempt to consider the weak magnetic field as 
a perturbation fails in the most interesting limit uj —>■ (one obtains terms like X 2 /uj). 
This reflects a clear physics: travelling a sufficiently long time the electron 'feels' any weak 
magnetic field. 
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The failure of the perturbation theory signals immediately that the effect of the weak 
magnetic field is more complicated than a slight changing of the eigenenergies of the effec- 
tive Schodinger equation and hence the localization length. The first excited state for the 
orthogonal ensemble has the eigenenergy S Q = 1 /2 and this energy definitely changes if the 
magnetic field is applied. But, as we will try to show below, an additional eigenstate with 
the eigenenergy £ u = 1/4, the energy of the unitary ensemble, appears (strictly speaking, 
a continuum of states with the energies £ p — (1 + p 2 )/4, where p is a continuous variable, 
should be added). Although at X <C 1 this state contributes to correlation functions with a 
small weight, the value £ u = 1/4 does not depend on X. 

As concerns the eigenenergy of the ground state, it is exactly zero due to the supersym- 
metry of the model. The eigenfunction \l/ of the ground state is a complicated function Q 
depending on the ratio X 2 juj. We did not manage to find it for an arbitrary value of this 
ratio and can write in the limits X 2 /uj — > (orthogonal ensemble) and X 2 /uj — > oo (the 
limit we are interested in now) only. Although the solution for the orthogonal ensemble is 
very well known in the standard parametrizationi'0, it has not been written before in the 
magnetic parametrization. To find it in the limit Aid, Ai c > 1 we use Eqs. (|36|, ^) taken 



at X = 0. The solution for the functions (i = 0, 1, 2) can be sought in the form 

*< (Aw, Ai c ) = \ ld Ri (z), z = \ ld \ic (52) 
Then, we obtain the following equation for Ri 

(H + Qz) R (z) = (53) 

{Ho + uz) Ri (z) + 2CozRq {z) = (54) 

(Ho + u>z) R 2 - ACjzRx (z) = (55) 

where 

«. = -^ (56) 

Eq. (|53|) is exactly the final equation obtained for the orthogonal ensemble in the standard 
parametrization. Writing the eigenfunction of the ground state \1/ in the form 

# = Aidi? (57) 

we find for R 

R = Rq (z (1 + 2 (7777* - ««*))) (58) 

In the opposite limit X 2 /u — > 00, the solution of Eqs. (|36| , j37j ) is completely different. 
Putting Cj = in this equations we find a solution 

^0 (Aic) = (1 + XX lc ) exp (-XAi c ) (59) 
= M> 2 = 
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Formally, one could write non-zero solutions for and ^2 proportional to \l/o- However, 
we understand that in the magnetic parametrization the ground state function \1/ should not 
contain at u = Grassmann variables because the effective Hamiltonian does not contain 
them. 

So, in the limit u — > the form of the ground state function \1/ changes discontinuously 
from Eqs. Q57|) to Eq. (|59|) as a magnetic field is applied. Physically, it may correspond to 
lifting degeneracies of some states of the initial electron problem by the magnetic field. 

The ground state solution, Eq. (|59|), is written for A^, Ai c 1. We have checked that 
this state survives and does not considerably change its form solving Eqs. ( |H3| ) in the limit 
uj = 0. The limit of the vanishing frequency considerably simplifies the equations because 
the "diffuson variable" Aid decouples from the "cooperon variables" Ai c and A c . Then, one 
can search for a solution \l/ depending on the cooperon variables only. We used the standard 
over-relaxation method with the Chebyshev acceleration to found numerically that up to 
X as large as 0.3, Eq. (|5TS|), describes the solution rather well. At higher values of X, no 
dramatic changes occurred to this state either. 

It is clear that in the limit X — > 00 essential values of the cooperon variables Ai c and A c 
are close to unity for an arbitrary uj and one obtains the function \l/ for the unitary ensemble. 
The solution found at u — can be used for non-zero frequencies provided the variable A^ 



is not very large. From Eqs. (36,37|) we can estimate the region of the applicability of the 
solution as 

Ai d Ai c cD < 1 (60) 

In order to find the density-density correlation function (x — x') and its moments 
Poo (x — x') at large distances we have to find not only the ground state but also at least 
the first excited state <f>\. After that one should calculate the proper "matrix elements" in 
Eq. Q23| ) or use Eq. fl2Ip. Although we argue that the eigenvalue £ u = 1/4 of the lowest 
excited state is exact, we are not able, as previously, to write the proper wave function in 
all regions of the variables exactly. The solution for X can be sought in the form, Eq. (p9|), 
corresponding to the symmetry of the block Q 12 of the supermatrix Q, which is the usual 
way of constructing the proper excited eigenstatesifl Explicit formulae can be written at 
J <C 1 provided the variables Aid, Ai c are large, A^, Xid 1, but limited from above by 
the inequality (|60|). In this limit one may keep the function fx only and we come to the 
following equation 

(nf - 2Ai^) fi P = £ (P) h P (61) 

Hx is given by Eq. (|37|) and p is a continuous real variable. The variables A^ and Ai c in 
Eq. fl61|) separate from each other and we write the solution f\ (see Eq. ( P9f) in the form 

/ip = ^ 1/2 Xp(lnA w )vI/ (A lc ) (62) 

with the function ^0, given by Eq. (|5^). The function x i v ) should be written as 

/2\ 1/2 

Xp (v)=l-\ sm[-(p/2)v + 6(p)} (63) 
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Eq. (|63D corresponds to a picture of "free motion" suggested in Ref.B for a description 
of the region of not very large Ai. In the case considered here, Eq. (p^), is applicable in 
the region determined by the inequality flSop . In calculation of asymptotics of the density- 
density correlator at large distances small values of p <C 1 give the main contribution and, 
in this limit, the phase 5 is also small, 5 (p) ~ p (see Ref.0). 
From Eq. (|6^) we obtain immediately 

£{p) = \(l + P 2 ) (64) 

The form of the solution fi, Eq. (|6"2|), differs from the ground state by a dependence on 
the diffuson variables A^. This is the same difference as the one for the unitary ensemble. 
It is not difficult to check that the variable Aid separates from the cooperon variables even 
if Aic is not large, Ai c ~ 1. Then, the solution fl62|) may still be used in the region specified 
Eq. ( |60D provided the function ty Q (Ai c ) is replaced by the solution ^ (depending only on 
the cooperon variables) of Eqs. (|32| , |33|) taken at uj = 0. 

Eq. (|62D with the function \1/ can be used also for an arbitrary X. In the limit X ^> 1, 
the function \1/ decays fast as the variables Ai c and A c deviate from 1. This means that the 
cooperon variables are 'frozen' and we come to the unitary ensemble. Then, the restriction 
(|60|) is not important and all formulae can be written for arbitrary AidcJ3t3. Of course, in 
order to prove the existence of the states with the eigenvalues S(p), (|62D, rigorously one 
should investigate the exact equations for arbitrary values of the cooperon and diffuson 
variables and prove that the solutions 'behave well' everywhere. It is a difficult task that 
can apparently be performed only numerically. 

Assuming nevertheless that this state exists and is given for Aid, Ai c ^> 1 by Eq. (|62|) 
(the inequality (|60D should also be fulfilled) we can estimate the proper matrix elements, Eq. 
([24]), determining the average density-density correlation function p^ (x) and its moments 
Poo (x), Eq. ([{]). However, trying to calculate the integrals in Eq. (^4j) one encounters 
a difficulty. The integrals over the matrix elements of Q can be reduced to integrals over 
the variables Aid, Ai c , A^. A c and the Grassmann variables. Proceeding in this way one 
should use the Berezinian J of the transformation to these variables, Eq. (p0|) , which is 
singular at Ai c = A c = 1. The singularity in the Berezinian leads usually to additional 
contributions. To avoid explicit calculations of these contributions it is more convenient to 
calculate not the physical quantities of interest themselves but their derivative over X. This 
leads to additional factors X\ c — X 2 C in integrands, thus compensating the singularities of the 
Berezinians. 

Calculation of the physical quantities using the spectral expansion, Eq. (^), is still very 
difficult. The alternative way of calculations is solving Eq. (fj^, and calculating the integral 
in Eq. (0). A simplified differential form of Eq. (^0|) can be written for the function 
determining the correlator p^ as 

W+^^W^oW (65) 

where k = kL cu . In Eq. (|65"D, we do not write Grassmann variables explicitly implying 
that is the non-compact central part similar to (1 — A^) 1 ^ 2 fi in Eq. (|29]). 

The derivative over the field (T nn ) x of the correlator T n , Eq. can be reduced to the 
form 
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dT nn jk) ^ ^d_ 
dX ~ dX 



J Xlf\^P^d\ ld d\ic (66) 



Writing Eq. (|66|) the Berezinian J, Eq. (0), was used in the limit Ai c , Aid 3> 1 and 

p(«) _ p( n ) _|_ p(«) 

Using the fact that the Hamiltonian Ti^) contains the variables Ay and Ai c separately 
let us expand the function (Ay, A lc ) in the eigenfunctions of the diffuson part of the 
Hamiltonian. Writing 

if } (Ay, = J2 X i^p On Ay) P£ } (AO (67) 
p 

we obtain the following equation for P$ (Ai c ) 

(£ (P) + <*) (Aic) + ^ 0) < } (Aic) = A (Aic) (68) 

where 



A (Ay) = 1 A" 7 3/2 Xp (In Ay) A? c * (Ai c ) ^A ld 



The operator is the part of the Hamiltonian Ti^ in Eq. fl3T|) acting on the variables 
Aic 

A similar expansion was done in Refill to consider the case of a strong disorder. The main 
contribution in the integral over Ay comes from large Ay and the integral formally diverges. 
When considering the unitary ensemble one had to cutoff the integral by Ay ~ 1/qj. At 



weak magnetic fields considered now, the cutoff can be determined from Eqs. ( |15| , |28|) . We 
see that the term with the frequency can be neglected if AidAi c u; <C 1, which leads to the 
cutoff A^ ~ (Ai c tl>) -1 . So, we estimate the function A (Ay) at p 1 as 

A (Aic) ~ (~tCo) 1/2 - n pXl^o (Ay) (69) 
In a similar way we can reduce Eq. (jB~6|) to the form 

dT nn (k) , ._si/ 2 _n d /" D (n)/ X s.j. /■> >. rfAi c 



- (-^) 1/2 - ^ / < } (Aic) *o (Aic) ^ (70) 

Now, we have to solve Eq. ( |68| ) and calculate the integral, Eq. ([701). In order to solve Eq. 
(IH) we introduce a Green function g of this equation 

(£ (p) + in) g kp (Ay, A' lc ) + Hfg kp (Ay, A' lc ) = A? c 5 (A lc - A' lc ) (71) 
and write the solution as 

P$ = [ 9k P (Aic, A'y) A (A' lc ) (72) 
Then, Eq. ( |70"D is reduced to the form 
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(/7 "" (/, * ) ' (-^r 2n ^ 1 9 kp (A lc) A'J * (A lc ) * (A'J ^ft 1 ^ ( 73 ) 



dX 



(AicA'ic 



The last step to be done is to find the function g and calculate the integral in Eq. (|73|). As 
usual,0@ the combination £ (p)-\-in enters the final formulae and one should perform Fourier 
transform to get the coordinate dependence. The point k = i£ (p) is a branching point and 
the integral over k can be shifted such that the integration is performed along the edges 
of the cut [i£ (p) ,ioo). The point i£ (0) = z/4 gives the value of the localization length 
whereas the integration over p leads to an additional power law prefactor |xp 3 ^ 2 . Being 
interested mainly in determining the exponential decay we can immediately understand 
that all correlators T nn decay as exp (— \x\ /4L CU ). To calculate the prefactor we may put 
n = i£ (p) in Eq. (|7T|) and find the proper solution g (Ai c , A' lc ), which satisfies the equation 

Hfg, (Xic, Kc) = Xfj (A lc - A' lc ) (74) 

and has to be substituted into Eq. (|73|). 

The solution of Eq. ([74]) can be found easily because one of the solutions \l/ (Ai c ), Eq- 
(pD|), of the corresponding homogeneous equation is known. We can write immediately the 
second solution $ (Ai c ) in the form 

$o (A lc ) = * (A lc ) dX (75) 

The function <3> (Ai c ) grows exponentially at Ai c — > oo and is proportional to Af c for small 
Ai c . Using the functions (Ai c ) and $ (Ai c ) we write the Green function g (Ai c , A' lc ) as 

n tx \ M _ J $ o (Ai c ) (A' lc ) , A lc < A' lc 

9o (A 1C) A lc ) - | ^ $o ^ A ^ ^ (76) 



Substituting the function g , Eqs. (|76| , |59| , |75|) , into Eq. (|73| ) we see that in the region 
1 <C Aic, A' lc < 1/1" the integral is logarithmic and is safely cut from below and above by 
the limits of this region. It is important that we remain in the region of large Ai c where our 
approximations are valid. As a result of the calculation we obtain for the moments p£} (x) 
of the density- density correlation function 



dp^ (x) 
dX 



Xln(l/X) |x|~ 3/2 exp(- \x\ /4L CU ) (77) 



It is interesting to note that the same prefactor determines correlations at coinciding points. 
Writing (0) as 



dpt ] (0) . 

.~\2n-i f (AidAic) n dXiddX 



dX J A ld A x 



lc 



V 2 

lc 



and integrating as previously over the region < A^ < (Ai c co>) 1 , 1 < Ai c < l/X we come 
to the result 

^f^~Xln(l/X) (78) 
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Eq. $n\ ) shows that the derivative over the magnetic field of all moments of the density- 
density correlation function decays with the localization length L cu of the unitary ensemble. 
The moments themselves can be obtained by integrating Eq. ([77]) over the magnetic field. 
The result can be written as 



■ 1 n 3/2 
/^'(■'■>* ("J 



exp I — ^— J + C n X 2 In ( — ) exp 



4L co J \XJ^\ 4L CU 



(79) 



where C n is a coefficient depending on n. The first term in Eq. flT9"|) is added assuming 
that at zero magnetic field one should have the standard result for the orthogonal ensemble. 
(In our previous worJll we wrote In 2 (1/X) instead of In (1/X) which was a result of a not 
sufficiently accurate estimate). The second term in Eq. ( [79"D becomes larger than the first 
one at distances exceeding 

x^~L c ln(l/X) (80) 

and decays in the same way for all the moments. As we have discussed in the previous 
Section, this should mean that this term is determined by rare splashes of the wave functions. 
These splashes can be due to a hybridization of states with eigenenergies very close to each 
other and it is not surprising that they are sensitive to the magnetic field. 

A more delicate question concerns the distribution of the logarithm of the wave functions, 
governing the behavior of "typical" wave functions. This quantity is necessary for comparison 
with a recent numerical study,0 where the logarithm of the transmittance was shown to 
behave smoothly between the orthogonal and unitary ensembles. 

Our new state with the eigenenergy £ = (1 + p 2 ) /4 contributes to the distribution of the 
logarithms also. Unfortunately, the procedure developed for calculation of the moments can- 
not be extended to a calculation of the distribution function in the most interesting relation 
of the variables t m exp (— \x\ / L cu ) which would give the distribution of the localization 
lengths. The problem is that one needs to know the analytical properties of the Green func- 
tion gk p , Eq. (|7lD , for arbitrary complex p. One can understand this on the example of the 
pure ensembles. In Eq. (|5l|) , computation of the distribution function required shifting the 
contour of the integration over p to the saddle point p = — i hxt/u ~ %. One could perform 
this procedure there because the integrand was known. 

In the crossover regime, one can rely only on the picture of the "free motion" applicable 
at small p, which make calculation of the distribution function not feasible. We can say 
only that the distribution function has an additional weight at t ~ exp (— \x \ /L cn ) so that 
the entire function is broader that for the pure ensembles. At the same time, the main 
body of the function is not known, which does not allow us to say much about "typical" 
wave functions. So, the results that we really have do not seem to contradict to the results 
of the numerical study of the average logarithms^, where the second length was not seen. 
They suggest also that the information the Borland conjecture refers to, the behavior of the 
averaged logarithm of the wave function, may be not sufficient to characterize the crossover. 
Instead, the entire distribution function is of interest since it can differ from the standard 
logarithmically normal even for x ^> L c . 

As concerns a numerical observation of the second length L cn in the moments of trans- 
mittances or density-density correlations, it may be very difficult because one needs to make 
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computation for a very large number of configurations to be able to make a reliable statis- 
tics of the rare events. In addition, a smearing due to the presence of leads can make the 
observation even more difficult. 



VII. DISCUSSION OF THE RESULTS 



We see from Eq. ([79]) that for weak magnetic fields the first term dominates at small 
distances. However, whatever weak the magnetic field is, the second term always prevails 
at sufficiently large distances. This is due to the fact that, although the small crossover 
parameter X enters the pre-exponential factor, the exponential function containing the 
localization length L cn = 2L co decays much slower. Comparing these two terms, we obtained 
the characteristic distance Xh, Eq- ( PD| ) where the both terms are of the same order of 
magnitude 

To explain Eq. (|7^) qualitatively one can follow argument suggested by Mott for de- 
scribing the AC conduct ivity0'H with minor modifications necessary in our case. 

Actually, the arguments presented below apply for localized states in any dimensionality 
and allow us to obtain the scale xh, Eq. fl8"0|), without any calculation. Naturally, the 
localization length L c in Eq. (|80|) depends on the dimensionality of the space. However, to 
be precise, we will speak in terms of the localization in wires. 

In quasi-one-dimensional samples all states are localized (in other words, J^°° p£> (x)dx = 
1). At the same time, for frequencies exceeding the energy E c = D /L 2 the system behaves 
like a metal. Therefore, eigenenergies of states localized in a region of the wire of size 
L c are separated by the mean level spacing of the order of E c . This picture allows us to 
estimate immediately the magnetic field H c at which the crossover between the orthogonal 
and unitary ensembles occurs. One should simply consider the region of the size L c as a 
closed grain and recall that the crossover between the orthogonal and unitary ensembles in 
a metallic grain occur ME at flux through the grain given by the relation 



E T (0/0o) 2 ~ A G , (81) 



where 0o is the flux quantum, A^r = (vV)' 1 is the mean level spacing in the grain (V is the 
volume), having the same order of magnitude as the Thouless energy Et = D /L 2 . For a 
quasi-one-dimensional grain of length L = L c one obtains the characteristic value X ~ 1, 
which means that the magnetic flux through the segment of the length L c is equal to the flux 
quantum O ^- At such fields the localization length of the main body of the wave functions 
changes considerably and this was clearly seen in the numerical work.0 

At the same time, we knowBi that the low frequency asymptotics of the level-level 
correlation function in an isolated grain changes from uj to uj 2 as soon as an arbitrarily weak 
magnetic field is applied. This corresponds to a finite probability of having two different 
levels at any small distance uj. These levels hybridize by the magnetic field and this leads 
to the low frequency asymptotics uj 2 typical for the unitary ensemble. 

The two-scale localization considered in the previous sections is of a similar origin, 
namely, one can always find two localized states that are arbitrarily close to each other 
in energy, although they may be separated by a large distance in space. Using this fact 
MottSS calculated the AC conductivity at finite frequencies uj. We do not study the con- 
ductivity itself but rather different kinetic quantities like e.g. the density-density correlation 
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function. In our case, the role of an external perturbation is taken over by the magnetic 
field. 

Following Mott's arguments ,011 let us try to find states whose energies slightly differ 
from the energy of a chosen state. As we have mentioned above, the typical separation in 
energy of states localized within a region of the length L c is of order E c . In neighboring 
regions of the lengths L c one cannot find such states either: Even if two states were occa- 
sionally closely located in two initially isolated systems, they would hybridize and split after 
bringing these two segments of the wire in contact. The value of the splitting would be E c 
again. However, if two states are separated by a large distance x 3> L c , the splitting energy 
A x decays exponentially 

A x ~E e exp(-x/L c ) (82) 

due to exponentially decaying overlap between the localized states. (Mott considered a 
tight-binding model; hence, the energy characterizing the splitting in his consideration was 
A x ~ Iq exp(— x/L c ), where Iq is of the order of band width.) 

Thus, one can find states a and b with an exponentially small energy difference consider- 
ing states localized sufficiently far from each other. Now, turning on an external perturbation 
(magnetic field in our case) we want to mix these states, which, like in isolated granules, 
would lead to a behavior typical for the unitary ensemble (with the doubled localization 
length L cu = 2L co ). However, in order to mix the states we need a sufficiently large matrix 
element A a f, of the vector potential between the states. If the states did not hybridize at 
all this matrix element would be exponentially small and this is similar the situation one 
encounters when calculating the AC conductivity at low frequencies. Mott0H suggested to 
take into account the states when they start hybridizing and we use this idea for our esti- 
mates. In this case the matrix element A a b is not exponentially small because the hybridized 
wave function is concentrated equally in the both localization centers. 

Using the modelling of the localization center in terms of an isolated granule of the length 
L c and applying the magnetic field H we obtain easily a characteristic energy window Eh 
of the order of 

E H ~ X 2 E C . (83) 

This estimate follows from Eq. ([X5|) by using the fact that in wires the energies E c , Et and 
Aq are of the same order. In an isolated grain states within this window are mixed and 
their correlations obey the statistics of the unitary ensemble. 
So, we come to the following picture: 

If two localized states a and b are separated by a distance x such that the corresponding 
overlap integral A x exceeds Eh, the states are not considerably affected by the magnetic 
field and its influence can be considered perturbatively. 

If two localized states a and b are separated by a very large distance x such that the 
overlap A x is much smaller than Eh, then the matrix A a ^ is exponentially small and the 
effect of the magnetic field on the states can be neglected. 

The magnetic field is important if the distance x between the localized states a and b is 
such that A x ~ Eh- For such states the matrix element A a b is not small but, at the same 
time, the levels are close to each other. The magnetic field H mixes the states a and b and 
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one can expect for the hybridized states the statistics of the unitary ensemble. Comparing 
the energies A x and E H 

exp(-x/L c ) ~ X 2 , (84) 

we obtain the characteristic distance Xu, Eq. (]80[) . 

Assuming that the mixing of the levels by the magnetic field leads to the correlations 
given by the unitary ensemble we come to the conclusion that the exponential decay of 
the density-density correlation function and other correlation functions characterizing wave 
functions should be determined at distances x > xh by the localization length L cu of the 
unitary ensemble. 

These simple arguments explain the two-scale behavior, Eq. (|79|), and give the correct 
second scale xh, although the doubling of the localization length can be obtained by the exact 
calculation only. At the same time, the probability to obtain states at a desired distance with 
eigenenergies very close to each other is small. This small probability may be compensated 
by a large contribution to moments of the density and therefore can be essential when 
calculating these quantities. As concerns calculation of the logarithm the large splashes are 
not crucial and their contribution is small due to the small probability, which is in agreement 
with the results of the numerical investigation of the transmittance£il(actually, the authors 
of this work suggested as one of possible explanation of the difference between the averaged 
logarithm and the moments existence of "anomalously" localized states). 

An additional suppression of the splashes can be a consequence of a smearing due to 
inelastic scattering and presence of leads. To observe the splashes one has to be at distances 
from the ends of the wire exceeding xu- At distances smaller that xu the weak magnetic 
field cannot have any considerable effect of the wave functions because the validity of the 
perturbation theory in H is determined by the value of the smearing and not by the energy 
separation of the relevant states. But the transmittance probes the wave functions at the 
end of the sample. In this region, the magnetic field can become important only if the 
relevant magnetic energy Eh becomes of the order of the smearing energy and the latter is 
of the order of E c near the ends. 

Recently, localization in disordered wires in a magnetic field was studied analyzing the 
autocorrelation function of spectral determinants (ASD)e9 and no indication of the two-scale 
localization has been found. In our opinion, this cannot be a great surprise because the ASD 
is a quantity that can give only a rough picture of what is going in reality. By definition, the 
ASD probes only correlations of the energy levels and not the spatial structure of the wave 
functions. Any localized states contribute to the ASD independent of the distance between 
the centers of localization. This makes hardly possible to say anything about the overlap of 
the localized states, which is so important for the existence of the second scale. We think 
that the fact that the second scale was not identified in Ref.0 is not in contradiction with 
our finding but rather demonstrates that the ASD is not a proper quantity to describe the 
effect involved. 

Let us now discuss how one can try to check experimentally our predictions. Considering 
the transmittance of single wire one should speak about the logarithm of the transmittance 
that is a self- averaging quantity. In this case, the rare splashes we discussed previously are 
not important and one should see only one localization length which is at weak magnetic 
fields close to the length L co of the orthogonal ensemble. However, one can measure the 
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conductance of a system of a large number of wires connected in parallel. This is exactly the 
set up of the experiment of Ref.00. The large number of wires in this case leads to averaging 
of the transmittance itself rather than of its logarithm. In such a situation one might see, in 
principle, the second scale L cu . Physically this means that although the splashes are rare, a 
wire, where such a splash occurs, gives the main contribution to the conductivity. 

Unfortunately, in the experimentllllizl, the localization length was not measured directly 
but was extracted from the dependence of the hopping conductivity on temperature. 

It is known that the conductance at very low temperatures is possible due to activations 
through Mott resonances, which leads to the picture of the variable range hopping. § The 
localization length L c enters directly the activation energy in this conductivity 



a = cr ex P 



To/T) i/(i +d ) t To „( vL f)-> . (85) 



where d is the dimensionality of the system. However, the localization length entering this 
formula is the length for a typical wave function and the splashes are not very important. 

As a matter of fact, Eq. (|85|) does not apply to ID samples. The conductivity of ID 
chains and wires should obey a simple Arrhenius-type instead 



a = <T exp (-T /T) , T ~ (vSL c ) . (86) 

Remarkably, even accounting for electron-electron interactions preserves this kind of 
dependence.^! It is this law that was discovered in recent experiments in the regime of 
strong localization00 and activation energy T was used to extract the localization length 
L c . Which localization length enters the conductivity in the ID case is not as clear as for 
higher dimensions and there could be a chance that the second scale can be seen here. Per- 
forming such experimental analysis at different temperatures one could try to observe the 
doubling of the localization length considered, because decreasing the temperature should 
result in a crossover from the activation behavior with the activation energy T in Eq. (jS6[), 
to the activation energy T /2. Indeed, at very low temperatures, T < (uSxh)~ , the electron 
hopping due to the overlap of far tails plays the dominant role, so that the larger localiza- 
tion length L cu should be used in Eq. fl8"E|). For temperatures {vSxh)~ < T < (uSLc)^ 1 , 
hopping between the main body of localized states is essential and the localization length 
L co should be substituted for L c in Eq. flgB|). 

Needless to say that the second scale may be observed only if the wires are sufficiently 
long such that the main contribution to the resistivity comes from the variable-range hopping 
inside the sample and the wire can be considered effectively as infinite. In this limit the 
resistivity should be proportional to the length of the wire. Measuring the conductance 
of a finite wire at ultralow temperatures where the resistivity grows exponentially with 
the sample length cannot help extracting the two-scale behavior as it has been discussed 
previously. 

VIII. CONCLUSIONS 

Studying localization properties of infinitely long disordered wires within supersymmetry 
technique combined with the transfer-matrix technique, we managed to calculate correlation 
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functions of interest. This procedure is based on working in the coordinate representation, 
Eq. (0), and generating any quantity of interest from the two-point correlation function, 
Eq. (0). Using this approach we obtained all moments and distribution functions of the 
density-density correlator and Landauer-like conductivity for the orthogonal and unitary 
ensembles, Eqs. (|(|f|§,^I). These expressions are exact and they are valid for arbitrary 
distances. 

Studying the crossover between the unitary and orthogonal ensembles for weak magnetic 
fields we have found moments of the density-density correlations at large distances x ^> L c . 
We have demonstrated that the far tail of these correlation functions characterizing the 
wave functions decays at arbitrary weak magnetic field with the localization length which is 
double as large as that of the main part. As all the moments decay with the same length, we 
argue that this behavior is due to rare splashes of the wave functions that are very sensitive 
to the magnetic field. In contrast to this, the localization length characterizing a typical 
wave function might exhibit a smooth crossover. 

We emphasize, however, that in the crossover regime, not only study of the typical 
localization length, but also of the whole distribution of the logarithm is of interest. It 
should differ from the standard log-normal one even for x > L c so long as x < xh, Eq. flS0|), 
being asymmetric and broad. At present, we can not predict how the entire distribution 
function of the averaged logarithm of the wave function will develop in the limit of infinitely 
long distances. 

Slightly revising the Mott arguments for the AC conductance, EE we have shown how 
one can extract the characteristic scale of the problem xh, without performing any calcula- 
tions. These arguments are quite general, which implies that our results may be valid for a 
disordered system of any dimension. 

From the exact treatment as well as from the qualitative arguments of the previous 
section we understand that the effect we have found is very sensitive to different kinds of 
the level smearing. Therefore, it may be observed at low temperatures only. 

As concerns a possibility of an experimental observation of the effect it might be seen at 
low temperatures in the regime of the hopping conductivity in a system of wires connected 
in parallel. The wires should be sufficiently long, such that the system would obey Ohm's 
law. Then, if the tails of the wave functions determine the hopping conductivity, one should 
expect the Arrhenius law but with different activation energies in the limit of very low 
temperatures. 

We hope that it can become possible to check our finding both numerically and experi- 
mentally, although this can be a difficult task due to a small probability of the wave functions 
with the high sensitivity to the magnetic field. 
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